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Abstract 



^f-s ' A noncolliding diffusion process is a conditional process of N independent 

CN ■ one-dimensional diffusion processes such that the particles never collide with 

each other. This process realizes an interacting particle system with long- 
p/ ' ranged strong repulsive forces acting between any pair of particles. When 

Ph ' the individual diffusion process is a one-dimensional Brownian motion, the 

r~| . noncolliding process is equivalent in distribution with the eigenvalue process 

jrt I oi an N X N Hermitian-matrix-valued process, which we call Dyson's model. 

For any deterministic initial configuration of N particles, distribution of par- 
ticle positions of the noncolliding Brownian motion on the real line at any 
— 1 , fixed time t > is a determinantal point process. We can prove that the 

^ ■ process is determinantal in the sense that the multi-time correlation function 

^^ . for any chosen series of times, which determines joint distributions at these 

(y-s I times, is also represented by a determinant. We study the asymptotic behav- 

^^ ' ior of the system, when the number of Brownian motions N in the system 

l/^ ■ tends to infinity. This problem is concerned with the random matrix theory 

^^ i on the asymptotics of eigenvalue distributions, when the matrix size becomes 

infinity. In the present paper, we introduce a variety of noncolliding diffusion 
processes by generalizing the noncolliding Brownian motion, some of which 
are temporally inhomogeneous. We report the results of our research project 
K> ' to construct and study finite and infinite particle systems with long-ranged 

j^ ' strong interactions realized by noncolliding processes. 
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1 Introduction 

In a system of A^ independent one-dimensional diffusion processes, if we impose a 
condition such that the particles never collide with each other, then we obtain an 
interacting particle system with long-ranged strong repulsive forces acting between 
any pair of particles. We call such a system a noncolliding diffusion process. 
In 1962 Dyson [13] showed that, when the individual diffusion process is a one- 
dimensional Brownian motion, the obtained noncolliding process, the noncollid- 
ing Brownian motion, is related to a matrix-valued process. He introduced a 
Hermitian-matrix- valued process having Brownian motions as its diagonal elements, 
and complex Brownian motions as off-diagonal elements. The size of the matrix is 
supposed to be A^ X A^. By virtue of the Hermitian property, all eigenvalues of the 
matrix are real, and Dyson derived a system of iV-simultaneous stochastic differential 
equations for the process of N eigenvalues on the real line R. In the present paper 
we call this stochastic process of eigenvalues Dyson's model. (Strictly speaking, 
it is a special case of Dyson's Brownian motion models with the parameter /3 = 2 
as explained below.) If we regard each eigenvalue as a particle position in one di- 
mension, Dyson's model is considered to be a one-dimensional system of interacting 
Brownian motions. Dyson showed that this system is nothing but the noncolliding 
Brownian motion [^ |23]. 

A probability distribution on the space of particle configurations is called a 
determinantal point process or a Fermion point process, if its correlation 
functions are generally represented by determinants [63 ESI EZ] • The noncolliding 
Brownian motion provides us examples of determinantal point processes: for any 
deterministic initial configuration of A^ particles, distribution of particle positions on 
R at any fixed time t > is a determinantal point process |15]. Moreover, by using 
the method developed by Eynard and Mehta for multi-layer random matrix models 
P^ I5U] , we can show that the multi-time correlation functions for any chosen series 
of times, which determine joint distributions at these times, are also represented by 
determinants [371 SSI \SM SS] • In the present paper we call such a stochastic process 
that any multi-time correlation function is given by a determinant a determinantal 
process [15] . 

We study the asymptotic behavior of the system, when the number of Brownian 
motions A^ in the system tends to infinity. Since, as explained above, the noncollid- 
ing Brownian motion can be realized by the eigenvalue process and the correlation 
functions are expressed by determinants of matrices, this problem is concerned with 
the asymptotics of eigenvalue distributions, when the matrix size becomes infinity. 
The latter problem is one of the main topics of the random matrix theory [SD]- In 
other words, our research project reported in this paper is to construct infinite par- 
ticle systems with long-ranged strong interactions by applying the results of recent 
development of the random matrix theory [371 SHI Sll HSl 1331 132] . 

In the present paper, we introduce a variety of noncolliding diffusion processes by 
generalizing the noncolliding Brownian motion. In Section 2 first we explain basic 



properties of diffusion processes treated in this paper, such as Brownian motions, 
Brownian bridges, absorbing Brownian motions, Bessel processes, Bessel bridges, 
and generahzed meanders. The transition probabihty density of a noncolhding dif- 
fusion process is expressed by a determinant of a matrix, each element of which is 
the transition probabihty density of the individual diffusion process in one dimen- 
sion (the Karlin-McGregor formula). This formula provides a useful tool for us to 
analyze noncolhding diffusion processes. In Section 3 we state the Karlin-McGregor 
formula and present basic properties of noncolhding diffusion processes. When such 
a Hermitian-matrix-valued process is given that its elements are one-dimensional 
diffusion processes, it will be a fundamental and interesting problem to determine a 
system of stochastic differential equations for eigenvalue process of the given matrix- 
valued process. Bru's theorem [HI IS] and its generalization [521 SD] give answers to 
this problem. In Section 4 we give a generalized version of Bru's theorem and show 
its applications. The determinantal structures of multi-time correlation functions of 
noncolhding processes are explained in Section 5. There asymptotics in A^ — )■ oo are 
also discussed [371 SH gSl El gS] . 

When we impose the noncolhding condition on a finite time-interval {0,T),T G 
(0, oo) instead of an infinite time-interval (0, oo), the noncolhding diffusion processes 
become temporally inhomogeneous, even if individual one-dimensional diffusion pro- 
cesses are temporally homogeneous. In Section 6 we discuss these temporally inho- 
mogeneous noncolhding processes. These processes are not determinantal any more, 
and make a new family of processes, which we call PfafRan processes [33, 112] ■ In 
the last section. Section 7, we list up the topics, which are related to noncolhding 
processes, but can not be discussed here. 

2 Brownian motion and its conditional processes 

Let {Q, J-", P) be a probability space. The stochastic process called (one-dimensional 
or linear) Brownian motion, {B{t,u)}t£[o^oo), satisfies the following conditions : 

1. B{0,u) = with probability one. 

2. For any fixed u E Q, B{t, u) is a real continuous function of t. (This property 
is expressed by saying that B{t) has a continuous path.) 

3. For any sequence of times, to = < ti < • • ■ < Im, M = 1,2,..., the in- 
crements {B{ti) — B(ti-i)}i=i^2,...,M are independent, and distribution of each 
increment is normal with mean zero and variance tj — tj_i. 

Then, the probability that the Brownian motion is observed in the interval [oj, bi] C 
R at time tj for each i = 1,2, . . . , M, P{B{ti) G [oj, bi\,i = 1,2, ... , M), is given by 



/ dxi j dx2---j dxMWG{ti-ti^i,Xi-Xi-i), 



where Xq = and 



x^ 



G{t, x) = ^-— exp ( -— ) , t > 0, X G R. 



The integral kernel G{s, x; t, y) = G{t—s, y—x) is called the transition probability 
density function of the Brownian motion. For any fixed s > 0, under the condition 
that B{s) is given, B{u),u < s and B{t),t > s are independent. This property is 
called a Markov property. A positive random variable r is called a Markov 
time, if the event {r < u} is determined by the behavior of the process until time 
u and independent of the behavior of the process after time u. The first time that 
a Brownian motion visits a given domain D, which is called the hitting time of D, 
is an example of a Markov time. In the definition of Markov property mentioned 
above, if a deterministic time s is replaced by a Markov time r, then we obtain 
a stronger property, called a strong Markov property. In general, a stochastic 
process, which has a strong Markov property and has a continuous path almost 
surely, is called a diffusion process. See [721 ES] for instance. In the case that the 
transition probability density function G{s, x; t, y) does not depend on times t and s 
themselves but only depends on the difference t — s, a. Markov process is said to be 
temporally homogeneous. In this case we write the transition probability density 
function as G{t — s,y\x) instead of G{s,x]t,y) to clarify its homogeneity in time 
in this paper. The Brownian motion is an example of a temporally homogeneous 
diffusion process. (It is also spatially homogeneous.) 

For d G N = {1,2,...}, using independent one- dimensional Brownian motions 
Bi{t),B2{t), . . . ,Bd{t),t > 0, a (i-dimensional Brownian motion is defined by a 
vector-valued diffusion process B(t) = (-Bi(t), B2(t), . . . , Bii(t)),t > 0. 

We want to consider the Brownian motion under the condition that it visits the 
origin at a given time T > 0. Since the probability that this condition is satisfied 
is zero, we first consider the Brownian motion under another condition such that it 
visits some point in an interval {—s, e) at time T, e > 0, and then define the original 
conditional process by taking the limit e I 0. The transition probability density 
function obtained in this limit is 

rTf , ^ G{T-t,0\y)G{t-s,y\x) . . ^, .^ ^^ 

G{s,x;t,y) = — — , < s < t < T, x,yeR. 

G{T — s,0\x) 

It is a temporally inhomogeneous diffusion process. We call this process a Brownian 
bridge of duration T and denote it by P'^{t), t G [0,T]. 

Although one-dimensional Brownian motion can visit any point of R, we consider 
the Brownian motion conditioned to stay positive forever. This conditional process 
Y(t), t G [0, oo) is temporally homogeneous process with transition probability 
density function G^^^'^\t,y\x); 

G^^/^\t,y\x) = -(G{t,y\x)-G{t,-y\x)\, t > 0, x>0,y>0, (2.1) 
G^'/'\t,y\0) = -yG{t,y\0), t > 0, y > 0. 



The distance of a three-dimensional Brownian motion from the origin, (i?i(t)^ + 
B2{tY + B3(tyy^'^,t > has exactly the same transition probability density as 
f l2.ip . and is called the three-dimensional Bessel process. In other words, the 
three-dimensional Bessel process Y{t),t G [0, oo) has two different representations, 
'the representation by a Brownian motion conditioned to stay positive' and 'the 
representation by a radial part of the three-dimensional Brownian motion'. We also 
note that Y{t) solves the following stochastic differential equation |62] . 

y(t) = Bit) + / —^ds, t > 0. 
Jo Y{s) 

Consider the Brownian motion X{t),t G [0,T] under the condition that it stays 
positive during a finite time-interval (0,T], with T G (0, oo). Then the conditional 
process is a temporally inhomogeneous diffusion process with the transition proba- 
bility density function Gj,' ' '(s,x;t,y); 



^(1/2,1) 



:s,x-t,y) = ^^rzY^{G{t - s,y\x) - G{t - s,-y\x)], (2.2) 



0<s<t<T, x>0, i/>0. 



G^^/'''\0,0;t,y) = ^^h{T~t,y)yG{t,y\0), tG (0,r], y > 0, 

where h{s, x), x > 0, s > Ois the probability that the Brownian motion starting from 
X > stays positive during the time-interval [0, s]. This conditional process is called 
a Brownian meander. Using two independent Brownian motions -Bi(t), B2{t) and 
a Brownian bridge /3^(t) of duration T, which is independent of Bi(t) and -82(^)5 we 
define a one-dimensional diffusion process by {Bi{ty + B2{t)'^ + fS'^ {t)'^y^'^ , t G [0,T]. 
We can prove that this process is identified with the Brownian meander. That is, 
the Brownian meander has also two different representations, 'the representation by 
a Brownian motion conditioned to stay positive during a finite time-interval (0,T]' 
and 'the representation by a radial part of the three-dimensional diffusion process 
(5i(t),52(t),/3^(t)),tG(0,T]'[71]. 

By comparing (12. ip with (12. 2p . we see that the distributions of the three-dimensional 
Bessel process Y{t) and the Brownian meander X{t), both starting from the origin, 
are absolutely continuous and satisfy 



PiX{.) G dw) = xT-^^fiXi-) e dw). (2.3) 



The equality (12. 3 p is called Imhof s relation 

The Brownian motion, which is killed at the origin, is called an absorbing 
Brownian motion in the domain (0, cxd). Let G{t — s, y\x) be the transition prob- 
ability density of this process. It is the density of the Brownian motion at time t, 
which starts from x > at time s (< t), restricted on the event that it stays positive 



in the time- interval [s,t]. The reflection principle of Brownian motion gives 

G{t — s, y\x) = G(t — s, y\x) — G{t — s, —y\x). 

The first formula in (12. ip means that the transformation of the transition prob- 
ability density G(t,y\x), given by {y/x)G(t,y\x), is identified with the transition 
probability density G^^^'^\t,y\x) of the three-dimensional Bessel process. It implies 
that the three-dimensional Bessel process Y{t) is the Doob h,-transforniation of 
the absorbing Brownian motion in the domain (0, oo). 

When d G N, the distance of the ci-dimensional Brownian motion from the origin, 
{Bi{t)'^ + B2{ty + ■ ■ ■ + Bd{tYY^'^, defines a one-dimensional diffusion process, which 
we call the d-dimensional Bessel process. The Bessel process can be extended to 
the cases with all positive real values of d as follows. With a parameter u G (—1, cxd), 
the transition probability density function of the 2(z/-|- l)-dimensional Bessel process 
Y^^^t), is given by 

G^-^\t,y\x) = ^7e^p(-^^)^^(T)' ^>0' ^>0'?/>0' 
G'^''\t,y\0) = —-/■ -— -exp -— , t > 0, y > 0, 

where T{z) is the Gamma function and Iu{z) is the modified Bessel function with 
parameter u [62]. The behavior of the Bessel process depends on the dimension d 
(the parameter v = [d — 2)/2). When d is greater than or equal to 2 (z/ > 0), the 
process has the origin as a transient point, and when d is less than 2 (—1 < z/ < 0), 
it has the origin as a recurrent point. Moreover, if and only if d is greater than or 
equal to 1 (i/ > —1/2), it is a semi-martingale [62] . 

Yor [71] introduced a family of diffusion processes with two parameters (i^, k), 
V G (— l,oo), K G (0,2(z/ -|- 1)), which includes the Brownian meander as a special 
case {v-iK) = (1/2,1), and he called each member of the family a generalized 
meander. The generalized meander X*^'^'''^(t), z/ G (— l,oo), kE (0,2(z/ + 1)), is the 
diffusion process with the transition probability density 

Gp^\s,x;t,y) = ^^l^^^G^^\t-s,y\x), 0<s<t<T, x>0,y>0, 
hrp' {s,x) 

Gp-\0At,y) = ^^^^^^^^ te{0,T], y>0, 

where hp^\t,x) = J^ dy G^'^^T - t,y\x)y-'', x > 0, t G (0,T] [H]. Then Imhof's 
relation (12.31) between the three-dimensional Bessel process Y(t) = Y^-^^'^\t) and the 
Brownian meander X{t) = X^^^'^'^\t) is generalized as 



P(X(-)(.) G dw) = ^/_[7^\9^ ^ 1 ^(^^^^(-^ ^ ^^) 
r(z^ + 1 — K/2) \ w{T) 



for the 2(z/ + l)-diinensional Bessel process Y^'^\t) and the generahzed meander 
X^'^''^\t). We remark that, though the parameter k of generahzed meanders is in 
(0, 2(z/+l)), we can discuss the cases k = and k, = 2(z/+l). The former corresponds 
to the Bessel processes, and the latter the Bessel bridges, which are the conditional 
Bessel processes to arrive at the origin at a fixed time T > 0. 

3 Noncolliding diffusion processes 

3.1 Karlin-McGregor formula 

In order to analyze noncolliding diffusion processes, it is useful to represent the tran- 
sition probability density functions by means of determinants. The representation 
is called the Karlin-McGregor formula [35] in probability theory, and the Lind- 
strom-Gessel-Viennot formula |1S1 EB EH] in combinatorics. It is also regarded 
as a stochastic-process version of the Slater determinant, which originally ex- 
presses a many-body wave function of free Fermion particles in quantum mechanics 



[Karlin-McGregor formula] ([Ml SHI EI]) Let G(s, x; t, y) be the transition 
probability density function of a one-dimensional diffusion process. On the line R 

set A^ starting points Xj, i = 1,2, ... ,N and A^ terminal points i/i, i = 1,2, ... ,N 
with xi < X2 < ■ ■ ■ < xn and yi < y2 < ■ ■ ■ < Vn, respectively. The transition 
probability density function of the system of A^ diffusion processes restricted on the 
event that they never collide with each other during the time-interval [s, t] is given 
by 

Go{s,x;t,y) = det (G{s,Xi;t,yj)). 

When N = 2, this claim is essentially equivalent to the reflection principle for a 
Brownian motion. That is, this formula can be regarded as a generalization of the 
reflection principle |41] . 

3.2 Noncolliding Brownian motions in a finite and an infi- 
nite time- intervals 

Let W;^ be a subset of R^ defined by W^ = {x E R^ : a^i < 0:2 < ■ • • < x^}, 
which is called the Weyl chamber of type A^r-i in representation theory [20] • The 
transition probability density of the absorbing Brownian motion in W^, that is, 
the density function of an A^- dimensional Brownian motion at time t, which starts 
from X G W^ at time 0, restricted on the event that it stays in W^ during the 
time-interval [0,t], is represented by 



fN{t,y\x)= det (G{t,yj\xi)) 

l<i,j<N \ J 



by the Karlin-McGregor formula. Then the probabihty that the Brownian motion 

tA 

N 



stays in W^ until time t is 



J\fNit,x)= / fNit,y\x)dy. 



Now we consider the noncolliding Brownian motion in a finite time- 
interval t e (0,T], X{t) = (Xi(t),X2(t),...Xjv(t)). The transition probability 
density of the process denoted by gN,T{s, x; t, y) is the conditional density of A^ Brow- 
nian motions at time t, which started from the points x = (xi,X2, . . . ,xi\f) G W^ 
at time s (< t), under the condition that they never collide with each other in the 
time- interval [s,t]. It is given by 



^.""ll ^''^\ fN{t-s,y\x), 0<s<t<T, x,yeW^. 



gN,T{s,x;t,y) = YTTT^ '—:fN{t-s,y\x), < s < t < T, x,y eWj. 

(3.1) 

The transition probability density in the case that all A^ particles start from the 
origin will be defined by taking the limit a:; — )• = (0, 0, . . . , 0) in fl3.ip . Since both 
of the numerator and the denominator in ( 13.1 p tend to as a; —^ 0, we have to know 
the asymptotic behavior of fN{t,y\x) and J\fN{t,x) in \x\/^/t — )■ 0. Performing 
bilinear expansions [38| HOj 113] with respect to the multivariate symmetric functions 
called the Schur functions l20l [T9l , we have obtained 



f~N{N+l)/4 / ^ \ / I |2\ 

fN{t,y\x) ~ c,(iv) ^^^TfJ^^^^^^^H'^J' ^^-^^ 

.r / ^ C2(N), f X\ \x\ 

Here hisf{x) is the N x N Vandermonde determinant, which is equal to the 
product of differences of variables Xi,X2, ■ ■ ■ ,xiy, 

hN{x) = ^ Jet ^ (x'-^J = Yl i^J - ^i)' 

and CiiN) = (27r)^/2n,1ir(^), C2iN) = 2^/^ n,1i r(z/2). By using ([X21) we ob- 
tain the transition probability density function of the noncolliding Brownian motion, 
when all A^ particles start from the origin {i.e. X{0) = 0) as 

rpN{N-l)/4^~-Ny2 / I |2\ 

gr^AO, 0; t, y) = ^-— Un{T - t, y)hM{y) exp ( "^ J > (3-3) 



tG(0,T], yew 



As illustrated in the left picture of Fig.l, in this case the A^ Brownian motions 
starting from the origin at time t = rapidly separate from each other to avoid 
collision. 
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Figure 1: Illustrations of the noncolliding Brownian motions X{t),t G (0, T] in the 
left picture, and Y{t),t G (0, cxd) in the right picture, both start from 

When the time T becomes infinity, the process X{t) converges to the temporally 
homogeneous process Y{t), the noncolliding Brownian motion in an infinite 

time- interval t G (0, oo), whose transition probability density function pi\i{t,y\x) 
is given by follows; 



PN{t,y\x) 



PNit,y\0) 



hNJy) 

hN{x) 



fN{t,y\x), t > 0, x,yeW 



A 

Af' 



(3.4) 



h^iyfexp 



\y\ 



t> 0, 



yew^. 



Ci{N) "'''"' """ V 2t 

The above formulas are derived form (13. ip and (13. 3 p by taking the limit T — ;■ oo 
using (13. 2p . See the right picture of Fig. 1, which illustrates Y(t) starting from 0. 
The first formula of (13. 4p implies that Y{t) is the Doob /i-transformation of the 
absorbing Brownian motion in W^, whose transition probability density is given 
by fN{t,y\x) [23]. One-parameter family of interacting Brownian motions on R 
satisfying the following system of stochastic differential equations 

2 



r,(t) = B,{t) 






Yd 



YA 



-ds, l<i< N 



(3.5) 



is called Dyson's Brownian motion model with parameter /3 > [HI |50]. It 
is readily seen from (13. 4 p that Y{t),t G (0, oo) solves the system of equations (13. 5p 
with (3 = 2. 

By comparing (12. 2p with (13.11) . (13. 3p . we find that Xit) can be regarded as a 
multi-dimensional extension of the Brownian meander X{t). Similarly, by comparing 
(12. ip with (13. 4p . Y{t) can be considered to be a mult i- dimensional version of the 
three-dimensional Bessel process Y{t). Moreover, a multi-dimensional extension of 
Imhof's relation ([O]) is derived from ([XI]), dSS]) and (El]) as [38j 

Ci(iV)T^(^-i)/4 



P{X{-) G dw) 



C2{N) h{w{T)) 



-P{Y{-) G dw). 



(3.6) 
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3.3 Noncolliding generalized meander and noncoUiding Bessel 
process 

We introduce the subsets W^ and W^ of R^ defined by 

W^ = {x e R^ : < Xi < a;2 < ■ ■ ■ < x^}, 
W^ = {a; G R^ : < |a;i| < ^2 < ■ ■ ■ < Xn}, 

wliicli are called the Wyle chambers of type Cn and of type Dn, respectively. The 
Bessel process Y^'^\t),t > has the origin as a transient point when z/ > 0, and it 
has the origin as a recurrent point when — 1 < z/ < 0. Then the state spaces Wjv of 
the noncolliding generalized meander and the noncolliding Bessel process, which will 
be introduced in this subsection, are W^ when u > 0, and W]^ when — 1 < z/ < 0. 
If the Bessel process is defined by the square root of the squared Bessel process, 
it can be a multi-valued stochastic process. Consider the squared Bessel process 
starting from a positive initial point. When z/ > 0, the process stays positive with 
probability one, and its square root is determined uniquely, which coincides with 
the Bessel process introduced in Section 2. While, when — 1 < z/ < 0, it hits the 
origin with probability one and then the square root process becomes a bi-valued 
process after hitting. The generalized meander and the leftmost particles in the 
N particle systems of the noncolliding generalized meander and of the noncolliding 
Bessel process are in the same situation. The absolute value |a:i| appearing in the 
definition of W]^ implies that bi-valued processes are allowed, when — 1 < z/ < 0. 
See Fig 2. However, we usually consider only nonnegative parts of such bi-valued 
processes just for simplicity of explanation. 

The density function of an iV-component generalized meander at time t, which 
starts from x in W^r at time s and stays in W^r up to time t, is given by 

fi:;^'\s, x; t, y) = det (g^^'''\s, x.; t, y,)) 

from the Karlin-McGregor formula. The probability that the process stays in W^v 
during the time-interval (0, t] is given by 



Afi,''^\t,x)= [ dyfi,''^\0,x;t,y). 

JWn 



Then the transition probability density function of the noncolliding generalized 

meander X('^''^)(t) = (X;"'"^(t),X^"'"^(t), . . .,xi;'"\t)) is given by 



AfS,''-\T-s,x 



g%p{s,x,t,y) = '-J-^^ '^fl;;'-\t~s,x,y), < s < t <T, x,y eW 



N- 



Consider the noncolliding generalized meander, when all A^ particles start from 
the origin. Define f}^ {t,y\x) = f^' '{0,x;t,y). Then we see 

/r«..w-r-n(^)exp{-ii:(x;+,?,l4et/'^^* 

10 





Figure 2: Illustrations of the noncolliding generalized meanders X'^'^'"^-' (t) with u > 
in the left picture, and — 1 < z/ < in the right picture. When — 1 < z/ < 0, a 
bi-valued process is assigned to describe the motion of the leftmost particle. 



A'^,'^) 



and using it g}^'j,'(s, x; t, y) is rewritten as 



A/jV '{T-s,x) 



g}^J{s,x;t,y) 



(3.7) 



where A/'jv^'''(t, x) = J^ dy f^ {t, y\x) Y\.i=i Hi '^- Again performing bilinear expan- 
sion with respect to the Schur functions, we have obtained the following asymptotics 
El: 



fVit,y\x) 



j.-N(N+l+2u)/2 / „ , 



(y) exp 



\y\ 



2t 



X 






^ /.i°W^l, ^^0, 



\x\ 



where C^'\N) = 2^(^+'^-i) [iti rWr(i + z/), C^''^\N) = 2^(^+2----i)/27r-^/2 
ni,{r(V2)r((z + 2z/ + 1 - fi:)/2)}, and 



AT 



4"(a) 



n ("I-"?)!! 

i<«<i<A^ 



af. 



fe=i 



By the above estimates the transition probability density function of the noncolliding 
generalized meander, when all A^ particles start from the origin is determined as 

rpN{N+K-l)/2A.-N{N+v) _ 

^i-)(0,O;t,y) = i —J- Ar(-)(T-t,y)/.r^^)(y)e-l?^l^n3.8) 



C^^^^){N) 



tG(o,r], yeWN. 
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Figure 3: Illustrations of the noncoUiding Bessel processes Y^'''(t) with z/ > in the 
left picture, and — 1 < z/ < in the right picture. When — 1 < z/ < 0, a bi-valued 
process is assigned to describe the motion of the leftmost particle. 

When the time T becomes infinity, the process ^^'^''^\t) converges to the tempo- 
rally homogeneous process Y^'^'(t), whose transition probability density function is 
given by 



p'}^\t,y\x) 



h 



(0) 



^fi:;\t,y\x), t>0, x,yeWN, 



j.-N(N+u) 

CH(Ar)^ l?/J exp 



(3.9) 



\y\ 



2t 



t > 0, y eW 



TV- 



See Fig 3. Since the parameter k controls the distribution of the process when 
t ^ T, it is irrelevant for the process Y^'^'{t) in which T — > oo. The process Y^'^\t) 
is the noncoUiding 2(z/ + l)-diniensional Bessel process, which is temporally 
homogeneous and solves the following system of stochastic differential equations 
when u > -1/2F 



Jo ^^[S) l<i<N-^0 V 



2Y}'\s) 






)2-yi"Vs)2 



ds, (3.10) 



t G (0,oo), 1 < i < A^, 



where we impose a reflecting wall at the origin when v = —1/2. Although solution of 
the system of equations is not necessarily unique in general, Y^'^\t) can be defined 
as a unique solution such that all coordinates are positive [&]. By comparing (13. 7p . 
(13. 8 p with (13. 9p . we have the following equality 

CMfp^) rpN(N + K-l)/2 

p(X('^'«)(.) G dw) = ^, ;, J, , , P(r('^)(-) G dw), 

C(-'«)(Ar) /iW(^(T)) 

which is an extension of Imhof's relation for the noncoUiding Bessel process Y^'^'(t) 
and the noncoUiding generalized meander X^'^''^'(t) ^0\ . 
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4 Matrix-valued processes 

4.1 Generalized Bru's theorem 

We denote the space oi N xN Herinitian matrices by 'H{N) and the spaces oi N x N 
real symmetric matrices by S{N). For a matrix A we indicate by *A the transposed 
matrix of A, by A the complex conjugate of A, and by A* = ^A the adjoint 
matrix of A, respectively. We denote the unit matrix of size A^ x A^ by I^. Bru 
[H IS] studied the eigenvalue processes of Wishart process, which is an 'H(A^)-valued 
process, and derived the stochastic differential equations for the eigenvalue processes. 
The result is generalized to the case that each element of matrix-valued process, 
^ijit), 1 < ^, J < ^, is a complex- valued continuous semi-martingale [39|H0]. In this 
section we state this generalized version of Bru's theorem and give its applications. 
Let A(t) = (Ai(t), A2(t), . . . , A7v(i)) be the vector, whose coordinates are eigenval- 
ues of 'H(A^)-valued process, E{t) = (^ij(t))i<i,j<Af, with Ai(t) < A2(t) < ■ ■ ■ < Aiv(i)- 
Then let U{t) = (Mij(t))i<i,j<Af be a unitary-matrix- valued process, which diagonal- 
izes S(t), 

U{tyE{t)U{t) = A(t) = diag(Ai(t), A2(t), . . . , A^(t)). 
We put 



T,,,k,{t)dt= (u{tydE{t)u{t)) iu{tydE{t)u{t) 

and the bounded variation part of {U{t)*d'E{t)U{t))ii is written as dTi{t). Then we 
introduce the Markov times 

aij = inf{t > : A,(t) 7^ A,(t)}, 

Tij = inf{t > aij : Xi(t) = Aj(t)}, r = min Tij. 

l<i<j<N 

[Generalized Bru's theorem] ([SIEIEHIIID]) Let ^y(t), 1 < i,j < N he 
complex-valued continuous semi-martingales. Then the eigenvalue process A(t) of 
S(t) solves the following system of stochastic differential equations: 

dXi{t) = dMiit) + dJiit), t G (0, r), l<i<N. 

where M{t) = (Mi(t), M2(t), . . . , M^(t)) is the martingale with dM,{t)dMj{t) = 
riijj{t)dt, and J(t) = (Ji(t), J2(t), . . . , JAr(t)) is the process with bounded variation 
given by 

dJi{t)= Yl Y7a ^7A'^{^^(t)^^At)}^^3Md^ + d'^i(^)■ 

i<j<iv^^W-Aj(tj 

Here 1{^} denotes an indicator function of a condition u. 
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Most of the examples shown in the next subsection are systems such that all 
particles start from the origin and rapidly separate from each other to avoid colli- 
sion. In these systems (Jjj = 0, 1 < i,j < N, and r = oo, that is, the repulsive 
forces among particles are strong enough to prevent any collision. For instance, in 
the stochastic differential equations of Dyson's Brownian motion models (13. 5p . the 
repulsive force becomes stronger as the parameter /3 becomes larger, and it is shown 
that r < oo, if < /3 < 1, and r = oo, if /3 > 1 [S3]. This corresponds to the 
fact that the Bessel process is transient, if the dimension d > 2 {u > 0), and it is 
recurrent, if < c? < 2 (-1 < u <0). 

4.2 Examples 

In this subsection we give examples of eigenvalue processes obtained by the gener- 
alized Bru's theorem. 

Let u e No = {0,1,2,...}, Bij{t), Bij{t), 1 < i < N + u,! < j < N he 
independent one-dimensional Brownian motions, and s{t) and a{t) be A^xA^ matrices 
whose element are given by 

if i < j, 
if i = j, 
if i > j. 



(i) GUE process 

Consider the 'H(A^)-valued process defined by S^^^(t) = s(t) + ^/^a{t), t G 
[0, oo). For any fixed t G [0, oo), S*^^^(t) is the 'H(A^)-valued random variable whose 
probability density function with respect to the volume element lA{dH) of Ti^N) is 

where TrA represents the trace of a matrix A, and Ci{N) = 2^/^7r^ /^. We denote 
the group oi N x N unitary matrices by U(A^). The probability fi'^^^{H, t)U{dH) is 
invariant under any unitary transformation H — )■ U*HU for any U G U(A^). In the 
random matrix theory, such a statistical ensemble of 'H(iV)- valued random variables 
is called the Gaussian unitary ensemble, GUE [50|l5T]. The probability density 
of eigenvalues of GUE is given by 





^i.«w, 


if i < j. 




^^B.M 


Sij{t) = < 


Bii{t), 
1 


if i = j. 


aij{t) = < 


0, 

1 ~ 


ctively. 


[ v^""^-^^'^' 


if i > j. 




[ '^'"^> 



t-^/2 , f xV ( \x\^ 



(7^^^(a3,t) = -— -/i^^ exp 



for X = {xi,X2, ■ ■ ■ ,xn) G W;^ [501 |5T]. Here Ci(A^) is the same constant as 
Ci{N) in (13.21) . By applying the generalized Bru's theorem to S^^^(t), we see that 
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A(t),t G (0,00) solves the system of stochastic differential equations of Dyson's 
Brownian motion model fl3.5p with /3 = 2. In Section 3.2 it was shown that the 
noncoUiding Brownian motion Y{t),t G (0, 00) solves the same equation. Then 
the equivalence in distribution of the noncoUiding Brownian motion Y{t) and the 
eigenvalue process A(t) of H*^^^(t) is established. 

(ii) GOE process 

Consider the iS(A^)-valued process defined by S*^°^(t) = s{t), t G [0, 00). For 
any fixed t G [0, 00), H^'-'^(t) is the iS(A^)-valued random variable whose probability 
density function with respect to the volume element V{dS) of S{N) is given by 

/i^°^(^,t) = — -^^^exp(^--Tr^2J, SeSiN), 

where C2{N) = 2^/^71^^^+^)/^. We denote the group of A^ x A^ real symmetric ma- 
trices by O(A^). The probability fi^'^^{S,t)V{dS) is invariant under any orthogonal 
transformation S -^^VSV for any V G O(A^). Such a statistical ensemble of S{N)- 
valued random variables is called the Gaussian orthogonal ensemble, GOE. 
The probability density of eigenvalues of GOE is given by 



^-N/2 



X \ \X 



|2 



for X = {xi,X2, ... ,Xj\f) G W;^ [SOI EI]- Here C2{N) is the same constant as 
C2{N) in fl3.2p . By applying the generalized Bru's theorem to S^^^(t), we see that 
\{t),t G (0, cxd) solves the system of stochastic differential equations of Dyson's 
Brownian motion model (13. 5p with /3 = 1. 

(iii) Laguerre process 

Let z/ G Nq. We denote the space of {N + z/) x iV complex matrices by Ai{N + 
v, N; C). Consider the M{N + u, N; C)-valued process defined by L{t) = {Bij{t) + 
y/—lBij{t))i<i<N+u,i<j<N. For any fixed t G [0, 00), L{t) is the M{N + u,N;C)- 
valued random variable whose density is given by 

f-N{N+u) /I \ 

/i^^^UE^^^ ^) ^ ——— exp ( --TrL*L j , L e MiN + u,N;C), 

where C3{N) = (27r)^*^^+'^). Such a statistical ensemble is called the chiral Gaus- 
sian unitary ensemble, chGUE [JT]. And the 'H(iV)-valued process defined by 
E:^(t) = L{t)*L{t), t G [0, 00) is called the Laguerre process |16]. The matrix E:^(t) 
is nonnegative definite and has only nonnegative eigenvalues. Applying the gener- 
alized Bru's theorem, we see that the eigenvalue process A(t) solves the following 
system of stochastic differential equations: 



Xr{t) 



Jo T^.^TArJo Ms) - ■^j[S) 



l<j<N • 
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and T = oo. Moreover, by Ito's formula we can prove that ^(t) = (fi:i(t), . . . , KAr(t)) = 
(A/Ai(t), . . . , y^XN(t)) solves the system of stochastic differential equations fl3.10p . 
which implies the equivalence of K,{t) with the noncolliding 2(z/ + l)-dimensional 
Bessel process Y^'^'(t). Since u is nonnegative integer for the chGUE and the La- 
guerre process, the dimension 2(z/ + 1) of the corresponding noncolliding Bessel 
process is positive and even ^6\ . 



(iv) Wishart process 

Let u G No and denote the space of (N+u) x A^ real matrices by Ai^N+u, N; R). 
Consider the A^(A^ + z/, A^; R)-valued process W(t) = (-Bjj(t))i<j<Af+i.,i<j<Ar. For 
any fixed t G [0, oo), W{t) is the Ai{N + u, N; R)-valued random variable, whose 
probability density function is given by 



chGOE 



.-N{N+u)/2 /I \ 

{W,t) = exp ( --Tr 'WW j , W G M{N + iy,N;R) 



where C4{N) = (27r)^^^"^'^)/^. Such a statistical ensemble is called the chiral Gaus- 
sian orthogonal ensemble, chGOE [51]. The iS(A^)-valued process defined by 
SW(t) = ^W{t)W{t), t G [0,oo) is called the Wishart process [9J. We can show 
that the eigenvalue process \(t) of the Wishart process S^(t) solves the following 
system of stochastic differential equations 

Jo ,^mJo \{s)-Xj{s) 






and r = oo. 



We are able to apply the generalized Bru's theorem to Hermitian-matrix-valued 
processes with additional symmetries. We introduce the matrix o"o and the Pauli 
spin matrices ai,i = 1, 2, 3 by 



1 o\ /on / -v/^n /I 



^° = '^ = >^0 i; ' ^^ = [l Oj ' ^^ = ^v^ )^'''= \0 -1 

Suppose that N > 2 and define 2N x 2N matrices Up = Ij^ ^ Cp, p = 0, 1, 2, 3. By 
definition Sq = hN- Let s''(t) = (sfj(t))i<ij<Ar, a''(t) = {a'^j{t))i<ij<N, < p < 3 be 
independent copies of s{t), a{t). By using them the 7{(2iV)-process S(t) is divided 
into 2x4 = 8 terms: 

2(t) = X^ I {s^{t) ® cTp) + (v^a^(t) ®ap)\. 

We consider the four Hermitian-matrix-valued processes represented by the following 
four terms, 

3 

^es{t) = J2{eee{t)®^p), = 1,2, e = ±, 

p=0 
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where 






laP{t), iie = l,p^3oi 9 = 2,p = 0, 
iie = l,p = 3oT9 = 2,p^0. 



Putting 'H0e{2N) = {H G "H i^Sg = £^eH}, we see that ^^^(t) takes values in 
'H0^{2N). Due to the symmetries of matrices, eigenvalues have the following prop- 
erties; 

(i) When e = +, they are pairwise degenerated; A = (wi, Ui, ijJ2, 002, • • • , oo^, ujn). 
(ii) When e = —, they are in the form A = (wi, —ui, 002, —0U2, • • • , ujn, —<j~'n)- 

(v) GSE process 

In the case of (0,6) = (2, +), a matrix S G 'H2+(2A^) is said to be a self-dual 
Hermitian matrix, if it has the symmetry *SS2 = S22 in addition to Hermitian 
property. The matrix can be diagonalized by a unitary-symplectic matrix. For any 
fixed t G (0,00) the statistical ensemble of S2+(t) is invariant under any unitary- 
symplectic transformation, and is called Gaussian symplectic ensemble, GSE. 
The probability density function of eigenvalues of GSE is given by (SHI EI] , 

where Cs{N) = (27r)^/^ ]^^^-^ r(22). The eigenvalues are pairwise degenerated and 
represented as A = {u!i,uji,uj2,^2, ■ ■ ■ ,^n,<^n)- This is known as the Kramers dou- 
blets in quantum mechanics. Applying the generalized Bru's theorem, we see that 
the distinct eigenvalues Wj, 1 < « < A^ solves the system of equations of Dyson's 
Brownian motion model (13.51) with /3 = 4, 

tu^{t) = B,{t) + y] /" ^ r^rfs, l<i<N. 

For a pair of degenerated eigenvalues X2i-i = A2J = oji , cr2i-i,2i = 00, 1 < i < N. 
All other pairs separately move and never coincide with each other, that is, r = 00. 

(vi) Matrix-valued process of class C 

In the case of {6, e) = (2, — ), a matrix S G 'H2-{2N) has the symmetry *S(t)S2 = 
— S2H(t) in addition to Hermitian property. We denote by sp(2A^; C) the Lie 
algebra of complex symplectic group represented by 2A^ x 2A^ matrices. Then 
'H2-{2N) ~ sp(2A^;C) H 1-L{2N). For fixed t G (0, 00) the statistical ensemble 
of H2- (t) coincides with the random matrix ensemble called class C introduced by 
Altland and Zirnbauer |[2j . The eigenvalues of a matrix in this class are in the form 
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A = {ui,—uji,uj2,—uj2,---,<^n,—<^n)- We denote the increasing sequence of non- 
negative eigenvalues by u!(t) = (uiit) , U2(t) , . . . ,C(;Ar(t)). By the generahzed Bru's 
theorem, we see that uj{t) solves the following system of equations 



UJi{t) 



Bi(t)+ / -^^ds+ Y^ / \^T-^ TT H ^-^ r^\ds, 

Jo '^i(s) ^^mJo l^iiV-^jis) UJi{s) + UJj{s) } 






\<i<N. 



It is also verified that r = oo and a;(t) G W^, ^t G (0, oo) with probability one. 
Comparing the above system of equations with fl3.10p . we can conclude the equiv- 
alence of a;(t) and the noncolliding three-dimensional Bessel process "K'^'^-'(t), in 
distribution |10]. Moreover, it coincides with the noncolliding Brownian motion 
under the condition that it never hits the wall at the origin |10] . 



(vii) Matrix-valued process of class D 

In the case of (^, e) = (1, — ), a matrix H G 'Hi-{2N) has the symmetry *H(t)Si = 
— SiS(t) in addition to Hermitian property. We denote by so(2A^; C) the complexi- 
fication of Lie algebra of special orthogonal group represented by 2N x 2N matrices. 
Then 'Hi-{2N) ~ so(2A^; C)n'H(2A^). For any fixed t G (0, oo) the statistical ensem- 
ble of Hi_(t) coincides with the random matrix ensemble called class D introduced 
by Altland and Zirnbauer [2]. The eigenvalues of a matrix in this class are also in 
the form A = (cji, — a;i,a;2, — W2, . . . , wat, —(jJn)- We denote the increasing sequence 
of nonnegative eigenvalues by uj{t) = {coiit) , ijj2(t) , . . . ,LL!iy(t)). By the generalized 
Bru's theorem, we see that uj{t) solves the following system of equations 

u,{t) = B,{t) + Y, / (^ rr + , .]_ , . ]ds, l<t<N. 

It is also verified that r = oo and u:{t) G W]^, ^t G (0, oo) with probability one. 
Comparing the above system of equations with (13.101) . we can conclude that uj(t) 
is equivalent in distribution with the noncolliding one-dimensional Bessel process 
Y^~^''^'{t) [10]. Since one-dimensional Bessel process is identified with a reflecting 
Brownian motion, oj{t) can be also regarded as the noncolliding reflecting Brownian 
motion HOl. 



5 Determinantal processes 

5.1 Fredholm determinant 

Let X be the space of countable subsets of R without accumulation points. For 
X = {xi,X2, ■ ■ ■ , xn) G R^ we write {x} for an element {xi, X2, ■ ■ ■ , xn} of X. For 
xn G R^ and A^' G {1, 2, . . . , A^}, we write Xn' for {xi, X2, ■ ■ ■ , xn') G R^ . For the 
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temporally homogeneous noncolliding Brownian motion Y{t), the A'- valued process 
^^{t) = {Y{t)} has the transition probability density function 

{PN(t - s,y\x), if s> 0, x,y e W^, 
PN{t, y\0), if s = 0, x = 0, ye W^, 

0, otherwise. 

We call the process C,'^{t) also the temporally homogeneous noncolliding Brownian 
motion in this paper. The A'-valued noncolliding Bessel process ^^''^(t) can be 
defined as well. For a sequence of time < ti < ■ ■ ■ < tj^ = ^ and a sequence of 
positive integers less than or equal to A^, {Nm}m=iJ the multi-time correlation 
function of C,'^{-) at {tm, {x}^J^}),m = 1,2, ... ,M, is given by 



Pn \ti, Xj^^; t2, Xj^^; 



I 



M TiN-N„ 



n^^.iR 



(M) 



M ^ N 



M-l 



J = A'm + l 



£=0 



,(0) 



where we put to = 0, x\^ = 0. Let Co(R) be the set of all real continuous functions 
with compact supports. For / = (/i, /s, . . . , /m) G Co(R)^, 6 = {9i,92,..., 9m) e 
R^, we put Xm{x) = e^""^""^"^^ - 1, 1 < m < M and x = {Xi,X2, ■ ■ ■ ,Xm)- The 
multi-time moment generating function 



^n{x;0) = e 



M 



N 



exp \^0mY^ fm{Xi^{tm)) 



.m=l im=l 



of C,^{t),t G [0, T] can be expanded by means of the multi-time correlation functions 
as follows: 



N N 



N M 



EE-En^Ln-;"Ln- 



Afi =0 A''2 =0 Nm =0 m= 1 

M N, 
X 

m=l i=l 






Nm 



R"M - 



Y\dx 



(M) 



i=l 



n Wxm[x^r^jpN [ti,x^^\;t2,x%l;...;tM,x''j^2) 



An A'-valued process ^(t) is called a determinantal process, if its multi-time 
moment generating function is written by Predholm determinant as 



^(X; 0) = Det Sm,nS{x -y) + K(t^, x; tn, y)Xn{y) 



(5.1) 



with a locally integrable function K. We call the function K the correlation kernel 
of the process. By definition of Fredholm determinant, the multi-time correlation 
function of ^{t) is then given by 



Pn I ti, Xj^_^, i2, x^^, . . . ,iMi Xj^^^ 



det 
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diffusions 


Wyle chambers 


matrix-valued pr. 


RM 


orth. poly. 


Brownian motion 


An -I 


GUE 

UiN) 


GUE 


Hn 


even-dim. Bessel pr. 


Cn 


Laguerre 


chGUE 


K, ^ G No 


3-dim. Bessel pr. 
(absorbing BM) 


Cn 


class C 

sp(2iV;C)n7^(2iV) 


class C 


1/2 


1-dim. Bessel pr. 
(reflecting BM) 


Dn 


class D 

so(2iV;C)n?^(2iV) 


class D 


-1/2 



Table 1: Noncolliding diffusion processes, random matrix ensembles, and orthogonal 
polynomials 



with a Y.m=i ^rn X Em=i ^m matrix 



Af^W ^(2) (M) 



■' J \<i 



<N^,l<j<N„l<rn,n<M 



The noncolliding Brownian motion ^^(t) is the determinantal process with the 
correlation kernel K.n'- 



N-l 



KNis,x;t,y) = < 



-y - 



n/2 



1 ^ ft 



2s ^\s 

n=N 



n/2 



X 



^r. 



2s 

X 



^n 



y 



/2s 



Vn 



if s < t, 



if s > t, 



where (fnix) = {y/n2"'n\} ^''^Hn{x)e ^ '^, n = 0, 1, 2, . . . , are the orthonormal func- 
tions on R associated with the Hermite polynomials Hn{x) |371ll3]. The noncolliding 
Bessel process ^^''^{t) is the determinantal process with the correlation kernel K]^ : 



^M 



K%'{s,x;t,y) = < 



n=0 ^ ^ ^ ' ^ 



E 

=0 

c 

E(; 



^y v- / M u (^\ V (y 



if s < t, 
if s > t, 



where (/?5;(x) = ^r(n + l)/r(i/ + n + l)x''/2L^(a;)e-^/^ n = 0, 1, 2, . . . , are the or- 
thonormal functions on R+ = {x G R : x > 0} associated with the Laguerre 
polynomials L^{x) with parameter v > —1 jl2]. See Table 1, which summarizes the 
correspondence between noncolliding diffusion processes with finite number of par- 
ticles, statistical ensembles of random matrices (RM), and orthogonal polynomials 
used to represent correlation kernels. 
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5.2 Scaling limits 

When the number of the diffusion process A^ goes to infinity, the asymptotic behav- 
iors of the noncolhding processes ^^ and ^^'^ are determined by the asymptotics of 
their correlation kernels K^r, K]^ in A^ — ;■ oo. Suppose that the correlation kernel 
converges under an appropriate scaling limit. Then the multi-time moment gener- 
ating function \I/Ar(x; 0) and the multi-time correlation functions p^ converge, and 
then the process converges in the sense of finite dimensional distributions. In the 
following, we discuss the bulk scaling limit and the soft-edge scaling limit for 
the noncolhding Brownian motion Y{t) and the hard-edge scaling limit for the 
noncolhding Bessel process Y^"'^ (t) [521 EQl EH El EZl [H SSI | 



1. [Bulk scaling limit] As A^ — )■ oo, the process ^'^{N + t) converges to 
the infinite-dimensional determinantal process, whose correlation kernel /C*^™ 
is expressed by using trigonometric functions. 



( 1 r^ 

- / rfMe(*-^)"'/2^os(M(x-y)), 

vr ./n 



r'''{s,x-t,y)={ 



sin(x — y) 



7i{x-y) 

r>00 



if s < t, 



if s = t, 



1 



due^'-'^^ /\os{uix-y)), ii s > t. 



2. [Soft-edge scaling limit] Define the scaled process Oa{N,t)^^{N^^^ + t) = 
{Fi(A^i/3 + t) -a(A^, t), Y2iN^/^ + t) - a{N, t), . . . , Yn{N^/^ + t) - a{N, t)} with 
a{N,t) = 2A^2/3 ^ jyi/3^ _ ^2^4_ As A^ -> oo, it converges to the infinite- 
dimensional determinantal process, whose correlation kernel A^^' is expressed 
by using the Airy function Ai(x), 



IC'''{s,x;t,y) 



du e(*-^)"/2Ai(a; _ u)Ai{y -u), if s < t, 

oo 

du e(*-')"/2Ai(x - u)ki{y -u), if s > t. 



3. [Hard-edge scaling limit] As A^ — )■ oo, ^^'"[N + t) converges to the 
infinite-dimensional determinantal process, whose correlation kernel /C '-''•* is 
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expressed by using the Bessel function Jy{x), 



K^'\s,x-t,y) 



du e(*-")"'/V^(ux)MJ^(u2/), 



2^{J,{2x)yJl{2y) - J,{2y)xJl{2x)} 

x"^ — y"^ 



if s < t, 



, if s = t, 



xy / du e^* ''^" ^'^J„{ux)uJ„{uy), if s > t. 



The above three infinite particle systems are all temporally homogeneous. The 
system obtained by the bulk scaling limit is spatially homogeneous, and the other 
systems obtained by the soft- and hard-edge scaling limits are spatially inhomoge- 
neous (see Table 2). These infinite particle systems are reversible and their equilib- 
rium measures are determinantal point processes [121 SS] • Osada [511 EZ] constructed 
diffusion processes whose equilibrium measures are determinantal point processes, 
by the Dirichlet form technique. Although the coincidence of Osada's processes and 
the above processes is expected |43j, it has not been proved yet. If the coincidence 
were proved, it would be concluded that the infinite particle systems obtained by 
the bulk scaling limit and the hard-edge scaling limit solve the stochastic differen- 
tial equations (13. 5p and (I3.10p with N = oo (see [58]). Nonequilibrium dynamics of 
determinantal processes with infinite numbers of particles have been studied, which 
show the relaxation processes to the stationary determinantal processes with the 
correlation kernels /C*^™ and /C^' [HI US] . There the theory of distributions of zeros 
and orders of growth of entire functions [17] is applied to analyze the determinantal 
structures of noncoUiding diffusion processes with infinite numbers of particles. 



5.3 Tracy- Widom distribution 

Consider the motion of the rightmost particle in the temporally homogeneous non- 
colliding Brownian motion C,'^{t) = {Y{t)}. For a fixed time t > 0, from (15. ip with 
M = l,ti = t,9i = 1, Xi{^) = '~l{a>a} the probability that the position of the 
rightmost particle is less than a G R is given by the Fredholm determinant as 



P max YAt) < a 



E 
Det 



N 



exp<^log(l{Y^(t)<«} 
6{x - y) - K7v(t, x; t, y)l{y>a} 



Then the distribution function F^g^^{a) of the position of the rightmost particle in 
the process obtained by the soft-edge scaling limit is given by 



i^max(a) = Det 6{x-y)-IC '{t,x;t,y)l{yya} 
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kernels (finite system) 


Hermite if„ 


Hermite H^ 


Laguerre L'^ 


kernels (infinite system) 


(bulk) 

trigonometric 

sin, cos 


(soft-edge ) 

Airy 

Ai 


(hard-edge) 
Bessel 


spatial homogeneity 


homogeneous 


inhomogeneous 


inhomogeneous 


Painleve equation 


Py 


Pu 


Pm 



Table 2: Determinantal processes and Painleve equations 



where the kernel /C^' was defined in the previous subsection. Tracy and Widom 
represented the distribution function (Tracy- Widom distribution) as 



{a 



exp 



[x — a)q[x 



)^dx 



with the solution q{x) of Painleve II (see for instance |54] j 



Pn 



(fq{x) 
dx'^ 



2q{x)^ + xq{x) 



satisfying the boundary condition q{x) ~ Ai(x), x — ?■ oo. For determinantal pro- 
cess, the distribution of position of the right-nearest particle to the fixed point (for 
instance the origin) is described by Fredholm determinant as well as the rightmost 
particle. In the bulk scaling limit [HT] and the hard-edge scaling limit [TIj the dis- 
tributions of right-nearest particles to fixed points are studied precisely, and are 
represented by solutions of Painleve V (Py) and Painleve III (-Pm), respectively (see 
Table 2). 



6 Temporally inhomogeneous processes 

The noncolliding processes discussed in the previous section are temporally homoge- 
neous diffusion processes, in which noncolliding conditions are imposed in the infinite 
time- intervals (0, oo). On the other hand, in Section 3 we explained that the non- 
colliding Brownian motion X{t),t G [0,T] is a temporally inhomogeneous diffusion 
process with transition probability density function (13. ip with (13. 3p . if noncolliding 
conditions are imposed during a finite time-interval (0,T]. Since J\fN{0,y) = 1, 
y e W^ by definition, the probability density function (I3.3p of the process X{T) 
coincides with g'^'^^ given by (14. ip . In other words, the distribution of the noncol- 
liding Brownian motion at the terminal time T of the noncolliding time-interval is 
equal to the eigenvalue distribution of GOE. When < t < T, however, the distri- 
bution of X{t) is different from the eigenvalue distribution of GOE. In particular, 
when < t <^ T, we see that it is close to the eigenvalue distribution of GUE by the 
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asymptotic behavior (13. 2 p of A/Af(t,y). /,From the above observations it is expected 
that X{t),t G [0,T] exhibits a transition from GUE to GOE as t approaches T. 
Remind that the off-diagonal elements of the GOE process S'~'*-*^(t) are one- 
dimensional Brownian motions and those of GUE process S*^^^(t) are complex 
Brownian motions. Hence, an 'H(A^)-valued process, whose eigenvalue process re- 
alizes X{t),t G (0,T], should satisfy the condition that each off-diagonal element 
behaves like a complex Brownian motion for < t <^ T, and it becomes to behave 
like a real Brownian motion as t /^ T [59j. We find that, if each imaginary part 
of off-diagonal element is given by the Brownian bridge of duration T, which was 
introduced in Section 2, this condition is fulfilled. Let Pjj{t), 1 < z < j < A^ be in- 
dependent Brownian bridges, which are assumed to be independent of the Brownian 
motions Bij[t), 1 < i < j < N used in the definition of Sij{t). Then we put 



afjit) 



0, if^=j 



and define an 'H(A^)- valued process by 



(t) = [, (t) + yZia5(t)) te[0,T]. [Q.l] 

V ■" / l<t.i<N 



Then we see that the eigenvalue process A (t) of E'^{t),t G (0, T] is a temporally in- 
homogeneous diffusion process and is equivalent in distribution with the temporally 
inhomogeneous noncolliding Brownian motion X{t),t G (0,T] with X{0) =0 [39] . 
The equivalence is proved by the fact that the eigenvalue process A^(t) of S-^(t) 
and that of S^^^(t) satisfy the generalized Imhof's relation f l3.6p . Remember that 
the eigenvalue process of S^^^(t) is identified with the temporally homogeneous 
noncolliding Brownian motion Y{t) in distribution as shown in Section 4. 

This result implies that the process X{t),t G (0,T] has two different represen- 
tations, 'the representation by a noncolliding Brownian motion with the transition 
probability density (13. 3p given by the Karlin-McGregor formula', and 'the repre- 
sentation by an eigenvalue process of S^(t) given by (16. ip '. This claim is a gen- 
eralization of the result that the three-dimensional Bessel process as well as the 
generalized meander have two different representations, 'the representation by con- 
ditional one- dimensional Brownian motions' and 'the representation by radial parts 
of three-dimensional diffusion processes'. 

The 'H(A^)-valued process 'EJ'it) is decomposed into an eigenvalue part A(t) and a 
unitary matrix part U{t). The latter representation of the process X'^{t),t G (0,T] 
implies that it is obtained from S-^(t) by integrating its unitary matrix part U{t). 
By this observation the following identity is derived [39], which is called the Harish- 
Chandra integral formula [21] or the Itzykson-Zuber integral formula |29j . 
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[Harish- Chandra integral formula] Let dU be the Haar measure of the 
space U(A^) normahzed as /uctv) ^^ ~ ^- ^'^^ ^ ~ {^iy^2, ■ ■ ■ , xjy) G W^ and y = 
(yi, ^2, • • • , Vn) e W^, put Ax = diag(xi, X2, . . . , xn) and Ky = diag(?/i, ^2, • • • , Vn)- 
Then for any cr G R the following identity holds: 

^jj^~TriA.-u*AyU)Vi2.-)^C^^N}a^ det (G{a',y,\x,)). 

U(JV) hNix)hN{y) ^<iJ<N \ '' ' } 

The above argument is also valid for the noncolliding generalized meander X.^^''^'{i). 
The matrix-valued process, whose elements are the complex-valued process having 
Brownian motions as its real part and Brownian bridges as its imaginary part, 



MT(t)=(i?,,(t) + v^/3;(t) 



l<j<Af+i/,l<j<Af' 



exhibits a transition from chGUE to chGOE. Its eigenvalue process is a complex- 
valued process and different from the noncolliding generalized meander. In stead of 
the process Mxit), we consider the 'H(iV)- valued process defined by 

SLW(t)=MT(t)*MT(t), te[0,T] 

The eigenvalue process A^W(t) = {Xf^it), Xl^{t), ..., X^^{t)) oiE }f{t)is the stochas- 
tic process with A^ nonnegative coordinates. We put KY^{t) = >/A^^(t), 1 < i < A^ 
and consider K'^'^{t) = {kY^ (t) , k^"^ (t) , . . . , k^ (t)) . Then the process ^^^(t) is 
temporally inhomogeneous and is identified with the noncolliding generalized mean- 
der X^''^''\t),t G (0,T] with z/ G No, K = z/ + 1 and X(''''')(0) = [40j. As another 
example, by setting (//,«;) = (1/2,1), we can construct an 'H(A^)-valued process, 
whose eigenvalue process is a noncolliding generalized meander exhibiting a transi- 
tion from class C to class CI. (For the definitions of the random matrix ensembles 
called class CI and class Dlll-odd/even mentioned below, see [751 El EOl ITO].) 

For a system of 2N independent Brownian motions, we impose the condition that 
pairs of (2z — l)-th and 2i-th particles meet at the terminal time T for 1 < i < A^, 
in addition to the noncolliding condition in the time-interval (0,T). Then we can 
show that the system realizes the eigenvalue process of the matrix-valued process, 
which exhibits a transition from CUE to GSE. 

Consider the noncolliding generalized meanders with {u, n) = (z/, i/ + l),i/ G 
No, {i^,n) = (1/2,0) and (z/, k) = (—1/2,0), with the above mentioned additional 
condition at t = T. We can prove that they realize the eigenvalue processes of the 
matrix-valued processes, which shows transitions from chGUE to chGSE, from 
D to class Dlll-odd and from class D to class Dlll-even, respectively |10] . 
(Here chGSE indicates the random matrix ensemble called the chiral Gaussian 
symplectic ensemble [51] ■) 

Several temporally inhomogeneous processes have two different representations, 
'the representation by noncolliding diffusion processes', and 'the representation by 
eigenvalue processes of matrix- value processes' (see Table 3). By identifying these 
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homogeneity 


homogeneous 


inhomogeneous 


1 dim. diffusion 


Brownian motion 


Brownian bridge 




Bessel process 


Bessel bridge 
generalized meander 


matrix- valued pr. 


GUE 


GUE-to-GOE, GUE-to-GSE 




chGUE 


chGUE-to-chGOE, chGUE-to-chGSE 




class C 


class C-to-class CI 




class D 


class D-to-class Dlll-odd 
class D-to-class Dlll-even 


process 


determinantal pr. 


Pfaffian process 


corr. func. 


determinant 


Pfaffian 


moment gen. func. 


Fredholm det. 


Fredholm Pfaffian 



Table 3: Temporally homogeneous and inhomogeneous processes 



two representations, Harish- Chandra (Itzykson-Zuber) formulas are derived for ma- 
trices with a variety of symmetries |10] . 

Recently, a family of stochastic processes, whose multi-time moment generating 
functions are represented by Fredholm Pfafiians [HIl 112] , has been intensively studied 
[67117]. We call such a stochastic process a PfafRan process. In general, the multi- 
time A^ point correlation function of a Pfaffian process is described by a Pfaffian of 
2N X 2N matrix (see Table 3). Since Pfaffians of 2N x 2N matrices are reduced 
to determinants oi N x N matrices in special cases, Pfaffian process is regarded as 
a generalization of determinantal process. Dyson introduced stochastic processes, 
whose A^ point correlation functions are represented by iV x A^ quaternion determi- 
nants, and they have been studied since then [151 UHl EI]- These processes are also 
members of Pfaffian processes, because quaternion determinants can be expressed 
by Pfaffians. We showed that the temporally inhomogeneous version of noncolliding 
Brownian motion X{t) and the noncolliding generalized meanders X^'^''^' are Pfaf- 
fian processes [371112]. By evaluating the asymptotics of Pfaffians in A^ — )■ oo, we can 
prove the existence of infinite-dimensional Pfaffian processes in appropriate scaling 
limits. They describe temporally inhomogeneous infinite particle systems. For the 
noncolliding generalized meander X^'^''^', the general form of correlation kernel is 
described by using Riemann-Liouville differintegrals 



7 Miscellanea 

1. In this paper, we have discussed noncolliding systems of one-dimensional dif- 
fusion processes in the unbounded domains, R and R+, which are related to Gaus- 
sian ensembles of random matrices. We can also consider noncolliding systems in 
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bounded domains. In particular, the systems on a circle have been studied and the 
relation with the statistical ensembles of random unitary matrices called circular 
ensembles are reported [131 [23 EHl ES]- For the systems on finite intervals, the 
transition probability density functions are described by using the Jacobi polyno- 
mials and the systems are related to the random matrix model called MANOVA 
(multivariate analysis of variance) model [12j. 

2. Dyson's Brownian motion models, which solve the system of equations (13.51) . 
form a family of processes with a parameter /3 > 0. In this paper by applying 
the generalized Bru's theorem we have clarified the correspondence between the 
eigenvalue processes associated with the random matrix ensembles GOE, GUE and 
GSE, and the Dyson's Brownian motion models with /3 = 1, 2 and 4. In particular, 
it was shown that, when P = 2, the process is also realized by the noncoUiding 
Brownian motion. Recently, a family of random matrix ensemble with a parameter 
/3 > is proposed, in which the eigenvalue distribution is give by 



|2 



where Cp{N) is the normalization constant. This random matrix ensemble is called 
the Gaussian beta ensemble, whose elements are tridiagonal matrices such that 
diagonal elements are independent Gaussian random variables and {k, /c+l)-elements 
and {k + 1, fc)-elements, 1 < k < N — 1, are independent random variables with %- 
square distribution with degree-of-freedom (A^ — k)(3 [12]. 



3. The ensembles of random matrices, whose elements are independent complex 
Gaussian random variables, is called the Ginibre ensemble [22]. The eigenvalues 
of matrices in this ensemble are complex in general and the probability density 
function is given by 

^^^°(-) = 7r4M n k-^.Pexp(-^), zee-, 

^Gin(^) ^^f^f^^ V ^ J 

where CGm(-^) is the normalization constant. Characterization of g'^^'^{z) has been 
intensively studied (see for instance [Mj 158]). 

4. In the present paper noncoUiding diffusion processes are discussed. NoncoUid- 
ing systems of discrete time Markov processes have been also studied. In particular, 
the system of noncoUiding random walks, called the vicious walk model [17], is an 
interesting and important model, since it is related to the representation theory of 
symmetry groups through the Young diagrams, the Young tableaux, and the Schur 

functions [SlElESlESlEaESllSSlElE]- 
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